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Abstract 

Charge state distributions of ®®Fe and ®®Ni projectile ions passing through thin 
carbon foils have been studied in the energy range of 1.44 - 2.69 MeV/u us¬ 
ing a novel method from the x-ray spectroscopy technique. Interestingly the 
charge state distribution in the bulk show Lorentzian behavior instead of usual 
Gaussian distribution. Further, different parameters of charge state distribution 
like mean charge state, distribution width and asymmetric parameter are deter¬ 
mined and compared with the empirical calculations and ETACHA predictions. 
It is found that the x-ray measurement technique is appropriate to determine 
the mean charge state right at the interaction zone or in the bulk. Interestingly, 
empirical formalism predicts much lower projectile mean charge states compare 
to x-ray measurements which clearly indicate multi-electron capture from the 
target surface. The ETACHA predictions and experimental results are found to 
be comparable for energies > 2 MeV/u. 

Keywords: X-ray spectroscopy, Electron capture and loss 

process, Charge state distribution. Ion-solid collision 


1. Introduction 

When an ion beam passes through the target a considerable amount of fluc¬ 
tuation takes place in charge state of the ion due to different atomic phenomena 
e.g. electron capture and loss processes, inner-shell ionization etc. occurring 
both in the bulk and at the target surface. In general, final charge state of 
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the ion depends not only on the parameters associated with incident ion (ini¬ 
tial ion velocity, charge and atomic number) but also on target characteristics 
(atomic number and density). Even though a monochromatic ion beam with a 
fixed charge state is passed through the medium, several charge states emerge 
out of the target However, after a large number of collisions, an equilib¬ 
rium in charge state distribution (CSD) is established, where certain balance 
in electron capture and loss processes are attained. The study of equilibrium 
conditions where charge state fractions (CSF) as well as mean charge state {q-m) 
reach to a certain stability are very crucial and have plenty of applications in 
various fields including atomic physics [5], nuclear physics [3], astrophysics [3], 
biophysics [S], energy loss experiments [5], accelerator designs [7], detectors [5]. 
Many experimental and theoretical groups have worked in these fields of research 
since 1950’s [9] and therefore many reviews or collection of data can be found in 
the literature based on electromagnetic methods and empirical formalisms (e.g. 
Allison [TU], Betz |TT], Wittkower [T3], and Shima [T1IT3]). 

Worth to mention that there are many experimental techniques like CRBS 
[T4] . recoil separator [El [16], TOF [E], electromagnetic methods [E] etc. which 
can be used to separate adjacent charge states and to study the charge state 
distribution of the respective atomic systems. The main demerit of these tech¬ 
niques is that they account for the total charge of the ion in the detectors placed 
a few meters away from the target. This implies that these techniques give a 
measure of electron-capture and -loss processes in bulk as well as at surface of 
the foil and cannot segregate the charge changing phenomenon occurring only 
in the bulk or of the surface. Therefore, these techniques are not appropriate 
to measure the charge state distribution right at the ion-solid interaction. This 
difficulty can be avoided by using charge less observables in the experiments to 
measure the charge states right at the ion-solid interaction. Interestingly, in the 
past several measurements have been carried out using x-ray spectroscopy to 
study various plasmas [Eiiinj. Nevertheless, the measurements with x-ray pho¬ 
ton detection have not yet been employed to study the CSD and other relevant 
parameters Qm, distribution width and asymmetric parameter of the projectile 
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ions during passage from any solid/gas targets. With this motivation, we confine 
the work to study the CSD and its parameters right at the ion-solid interaction 
zone using the x-ray spectroscopy technique. 

2. Experiments 

The experiments were performed with the energetic ion beams of ®®Fe and 
®®Ni using 15 UD Pelletron [21] accelerator at lUAC, New Delhi. Well-collimated 
ion beam in the energy range of 1.44 - 2.69 MeV/u were bombarded on 80 
/ig/cm^ (« 113 /ig/cm^ at 45°) thick amorphous carbon target foils to produce 
the equilibrium charge state distribution. The target was placed at 45° to the 
beam axis so that the x-ray spectra could be measured right from the ion-solid 
interaction zone. The x-rays were detected in a Low Energy Germanium Detec¬ 
tor (GUL0035, Ganberra Inc., with 25 ^m thick Be entrance window, resolution 
150 eV at 5.9 keV) placed at 90° to the beam axis to avoid the Doppler shift. 
The x-ray produced in the ion-solid interactions were passed through two col¬ 
limators of 3 mm diameter kept at 55 cm apart whereas the first collimator 
was placed at 10 cm away from the target. This configuration ensured that the 
x-rays were coming from a tiny section (±3-4 mm) of the ion-solid interaction 
zone. In the time scale, the x-ray detector could observe only atomic transi¬ 
tions of very short life-time (few tens of psec) with respect to the centre of the 
interaction zone. Hence, the x-ray spectroscopy technique could be considered 
as a measurement right at the ion-solid interaction zone compared to the elec¬ 
tromagnetic measurements taking place away from the interaction zone or at t 
« a few /isec for MeV ions. The x-ray detector was placed outside the chamber 
at 65 cm away from the target separating a thin mylar window of 6 ^m at the 
interface of detector and chamber. The beam was dumped in a Faraday cup. 
Two silicon surface barrier detectors were used at ±10° to monitor the beam 
direction. The x-ray spectra observed for all the beam energies are shown in 
Fig. 1. Galibrations were done for the x-ray detectors using ®°Go and ^^^Am 
standard radioactive sources. The resolution was found to be about 200 eV 
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at 6.41 keV with the experimental conditions in the beam hall. Further, the 
calibration was internally verified through Fe and Kp peaks due to beam 
halo hitting the carbon foil holder made up of stainless steel in the case of ®®Ni 
projectile ions. However, in the case of ^®Fe-beam experiment, beam halo was 
minimized by passing the beam through a blank target frame so that its pres¬ 
ence did not affect much the peak structure originated from the projectile ions. 
Vacuum chamber was maintained at a pressure around 1x10“® Torr. 



Figure 1: X-ray spectra for (a) ®®Ni beam and (b) ®®Fe beam on 80 ^g/cm^ C-foil at different 
beam energies and initial charge states 


3. Data Analysis and Results 

In this work we are intended to determine the charge state distribntion of 
projectile ions right at the ion-solid interaction along with and other relevant 
parameters from the measured x-ray spectra. Accordingly we develop a novel 
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method to extract required information from the x-ray spectra observed. Worth 
mentioning here that the parameters like charge state fraction, mean charge 
state etc. obtained at different energies for the particular ion can be compared 
without normalizing the x-ray spectra. Hence, normalization of x-ray spectra is 
not carried out in this work like any other electromagnetic methods coupled with 
position sensitive detectors [22]. It is clear from the spectrum of 136 MeV ®®Ni 



Energy (keV) Energy (keV) 


Figure 2: X-ray spectrum of ®®Ni on C at 136 MeV (a) Fitting shows only broad features of 
the spectrum (b) Projectile x-ray peak is fitted into nine Gaussian functions corresponding to 
x-ray lines appearing from H-like to F-like ®®Ni with a exponential function representing the 
background (c) The residuals of fitting (b) is shown 

on C that it contains mainly three structures as shown in Fig. [^and]^. In our 
earlier work the first peak (7 - 8.4 keV) is recognized to have originated from the 
projectile ion x-ray, whereas second and third peak (8.4 - 11.5 keV) belong to the 
projectile-like fragment ions emanating from the nuclear reactions, respectively 
Ham]. It is worth mentioning here that we are only considering the charge 
changing phenomena in elastic events or in the projectile ions, thus the second 
and third structures are of no relevance in this work; hence they will not be 
brought in the further discussion. Another point is to be cleared at this stage 
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that when the fast projectile ions incident on target atom, the collisions create 
vacancies in the different shells of both target atoms and projectile ions, which 
cause shift in the characteristic x-ray line energies 1211 • static target 




Figure 3: Comparison between experimental and ETACHA charge state distribution for (a) 
^®Ni on C (b) ^^Fe on C for various beam energies. Figure is showing Lorentzian fit (red) to 
experimental data and Gaussian fitting (blue) to ETACHA data [271). Due to small values, 
errors are embedded in symbol itself. 


atoms are completely filled before the collisions, so the shift in K-shell x-ray 
energies depends on vacancies created in K-, L-, M- shells etc. during ion-solid 
interactions. In contrast, in the case of projectile ions the outer electrons get 
stripped completely up to a certain level during ion-solid collision (depending 
upon the collision energy) [5^ and the line energy shift solely depends on K- 
and L-shell vacancy only due to absence of electrons in the M-, N- and further 
shells [25]. The increased ionic charge of the projectile ions shift the centroid 
towards higher energy side with the increase in beam energies. Such a scenario 
is clearly depicted in Fig. 
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Important to note that life-time of a single K-vacancy is of the order of few 
fsec. Whereas in the intermediate energy regime (this work), projectile ion takes 
a few tens of fsec before leaving the target thus encountering multiple collisions. 
It causes a complex chain of collision events resulting in a composition of many 
charge states under the projectile x-ray peak. Though poor detector resolution 
constrains us to resolving the individual x-ray lines but the well-defined centroid 
gives a measure of the mean charge state of the projectile ions. To find the mean 
charge state and charge state distribution; projectile x-ray peak has been fitted 
with various Gaussian peaks ranging centroid energy for H-like to F-like ions, 
as shown in Fig. 2(b). It should be noted here that x-ray energies for Ni are 
available only for H-like to Li-like Ni [55]. The rest of the energies for Be-like 
to F-like Ni are scaled from corresponding Fe data [28] . 

In this way we have obtained the distribution of CSFs (Fg) due to K-shell 
vacancy production directly from the measured x-ray line intensities as follows 


^ IgCUg 

^ y I 


q J-q^q 


( 1 ) 


here, Ig and ujg corresponds to intensity and fluorescence yield of the projectile 
ion in the final electronic configuration, respectively. Noteworthy that, for the 
values of fluorescence yield, one needs to know the vacancy distribution in the 
initial configuration state, corresponding sets of the transitions schemes and 
the transition rates during heavy ion collisions. This can be easily estimated 
for target atoms. However, in the case of projectile ions (this case) neither is 
very well known nor is an easy problem to solve. Recently, Dr. M.F. Hasogulu 
has theoretically calculated LSJ-state dependent fluorescence yields in the case 
of Li-like to F-like ions in the range of 3<Z<30 atomic numbers for single K- 
vacancy cases [55]. We have taken the fluorescence yield data for different charge 
states from their calculations to estimate the fluorescence yield for a particular 
ionic state by statistical averaging of most probable LSJ states during ion-solid 
collisions, as discussed in detail elsewhere [SD]. The errors given for the CSF data 
is a combination of the systematical error (estimated, 5%) and the statistical 
error of the measured CSF. This procedure is repeated for all beam energies for 
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both the ion species. All the measured as well as ETACHA predicted charge 
state distributions at each beam energy have been plotted in Fig. Worthy to 
mention that the CSDs predicted by ETACHA m follow the Gaussian shape, 
in contrast the experimental values depict a different pattern, which is fitted well 
with a Lorentzian function for all the cases. Observation of Lorentzian profile 
for the charge state distribution is an important characteristic of any plasma 
[SI ED- This suggests that ion-solid collisions constitute a tenuous high density 
plasma in the bulk of the solid target, which is described in detail elsewhere 

m- 

Having known the charge state distribution, we step forward to measure the 
charge state distribution parameters also. Conventionally, the qm is computed 
using the CSF distribution as follows [31] 

Qm = Y^ qFq ( 2 ) 
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nevertheless, we have adopted a different approach here. Since the centroid 
of the Lorentzian distribution represents qm, we have obtained corresponding 
gm’s from the fitting as shown in Fig. [^and their uncertainties are simply the 
fitting error of each centroid, Aq^- Measured are compared to ETACHA 
calculations m shown in Fig. It is quite clear from the figure that ®®Ni data 
show some departure from ETACHA predictions m, however they agree quite 
well with each other in case of ^®Fe. Next, for comparison with the available 
empirical formalisms, we have choose Schiwietz formalism [35] because of its 
vast and updated dataset and compared it with experimental results, shown in 
the same figure. It is found that the empirical formalism shows lower values 
than both experimental values and ETACHA predictions. 

This is due to the fact that the highly charged ions with the low and interme¬ 
diate incident energies interact strongly with the target surface. In the earlier 
work it is reported that projectile ions captures target surface electrons into 
high Rydberg states to form hollow atoms while entering the target [351 EZl EH] • 
However, due to the low energy only the interaction with the entrance surface 
could be studied so far [39] . In the intermediate energy region (this work) when 


the projectile ions leaves the target surface, they can again captures surface elec¬ 
trons into high Rydberg states. However, due to long decay-time (of the order of 
/xsec-nsec m) these transitions could not be observed using x-ray spectroscopy 
technique as discussed earlier. Whereas, traditional electromagnetic techniques 
which are involved in measuring CSD and q^, are mainly electromagnetic in 
origin, thereby accounting for the total charge of the ion in the detectors placed 
at the focal plane, a few meters away from the target. Empirical formalisms 
are tuned on the basis of the measured data taken from the electromagnetic 
methods and therefore represent the integral role of the bulk and the surface 
of the foil. This is the reason why there is a long discrepancy occurs between 
results of empirical formalism and x-ray method in the low and intermediate 
energy regime, as reported earlier also mnn]. 

Whereas, ETACHA formalism m takes account of ionization and capture pro¬ 
cesses theoretically and does not include surface electron capture processes. 
Therefore ETACHA predictions ought to be represent the measurements right 
at the ion-solid interaction and may compare well with the experimental results 
(this work, x-ray method) in the intermediate energy regime. It is worth to note 
that as the projectile energy increases probability of capturing electrons from 
the target surface start decreasing and at much higher energies (GeV) due to 
very low interaction time, projectile will not capture electrons from the target 
surface. Thus at higher energies (GANIL energies) x-ray method, ETACHA 
predictions and empirical formalism or electromagnetic methods will start fol¬ 
lowing same results due to exclusion of surface effect. 

Next we have evaluated distribution width from the given formula [34] 

dP = Y,{q-qmfF, (3) 

and plotted them against beam energies in Fig. Uncertainties in distribution 
width are determined using standard procedure of propagation of errors. The 
figure shows good agreement between the measured and ETACHA predicted 
distribution widths for ®®Ni whereas they depart from each other with higher 
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Figure 4: Comparison between experimental, ETACHA and Schiwietz mean charge state for 
(a) ®®Ni on C and (b) ®®Fe on C on 80 /rg/cm^ C-foil. Solid lines are to guide eye only. Error 
bars are smeared within the symbol size. 

beam energies for ®®Fe. However, the oscillatory nature in distribution widths 
is quite common as observed earlier, for example [41) . 




Figure 5: Comparison between experimental and ETACHA distribution widths for (a) ®®Ni 
on C and (b) ®®Fe on C 

Further, we have computed the asymmetry of the charge distribution using 
the formula given by |34j 

s = ^iq-qmfFg/d^ (4) 

and the data so obtained for experiments as well as ETACHA predictions are 
plotted in Fig. Again, using standard propagation of errors procedure uncer¬ 
tainties are estimated in the asymmetry of the charge distribution parameters. 
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Figure 6: Comparison between experimental and ETACHA asymmetry of charge distribution 
parameter for (a) ^^Ni on C and (b) ^®Fe on C 

In both the cases of ®®Ni and ^®Fe a decreasing trend with increasing energy is 
found similar to the ETACHA values. 

4. Conclusion 

In this work we have determined the charge state distribution and its dif¬ 
ferent parameters like mean charge state, distribution width and asymmetry 
parameter using x-ray spectroscopy technique. This technique is found to be 
appropriate to segregate the charge state distribution in the bulk from that of 
the surface by measuring the charge changing phenomena right at the interac¬ 
tion zone. It is shown that ETACHA code m represents well the mean charge 
state measurements starting from energies > 2 MeV/u. However, an unusual 
charge state distribution is observed in the form of Lorentzian distribution in 
contrast to the prevailing Gaussian distribution as predicted by the ETACHA 
code. Whereas other charge state distribution parameters are found to show¬ 
ing comparable results to ETACHA predictions [27]. It is observed that the 
Qra data from either this measurement or ETACHA predictions m are much 
higher than that from any empirical formula, which is a clear indication of the 
multi-electron capture from target surface as reported earlier (HU 03103] . 


II 
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